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Abstract 

We report on resuhs from Suzaku broadband X-ray observations of the southwest part of the Galactic 
supernova remnant (SNR) RX J1713.7— 3946 with an energy coverage of 0.4-40 keV. The X-ray spectrum, 
presumably of synchrotron origin, is known to be completely lineless, making this SNR ideally suited for 
a detailed study of the X-ray spectral shape formed through efficient particle acceleration at high speed 
shocks. With a sensitive hard X-ray measurement from the HXD PIN on board Suzaku, we determine the 
hard X-ray spectrum in the 12-40 keV range to be described by a power law with photon index P = 3.2 ± 0.2, 
significantly steeper than the soft X-ray index of P = 2.4 ± 0.05 measured previously with ASCA and other 
missions. We find that a simple power law fails to describe the full spectral range of 0.4-40 keV and 
instead a power-law with an exponential cutoff with hard index P = 1.50 ±0.09 and high-energy cutoff 
€c ~ 1.2 ±0.3 keV formally provides an excellent fit over the full bandpass. If we use the so-called SRCUT 
model, as an alternative model, it gives the best-fit roUoff energy of eioii = 0.95 ±0.04 keV. Together with 
the TeV 7-ray spectrum ranging from 0.3 to 100 TeV obtained recently by HESS observations, our Suzaku 
observations of RX J1713.7— 3946 provide stringent constraints on the highest energy particles accelerated 
in a supernova shock. 

Key words: acceleration of particles — ISM: individual(RXJ 1713.7—3946) — ISM: supernova rem- 
nants — X-rays: ISM 



1. Introduction 

Over the past decade, the spectral and spatial prop- 
erties of non-thermal X-ray emission, presumably syn- 
chrotron radiation by multi-TeV electrons, in young shell- 
type supernova remnants (SNRs) have been extensively 
studied with X-ray satellites, especially with ASCA and 
Chandra, providing observational support for the gen- 
eral picture that Galactic cosmic rays arc sustained by 
strong shocks in SNRs (Koyama et al. 1995; Slane et al. 
1999; Gotthelf et al. 2001; Hwang et al. 2002; Uchiyama 
et al. 2003; Long et al. 2003; Bamba et al. 2003; Warren 
et al. 2005). An excellent example comes from Chandra 



observations of the Tycho's SNR, where the remnant is 
delineated by a thin (1" or less) layer of X-ray syn- 
chrotron emission that exhibits spectral steepening be- 
hind the shock (Cassam-Chenai' et al. 2007) . This strongly 
indicates that TeV-scale electrons are indeed accelerated 
at the outer blast wave. The well developed theory of 
diffusive shock acceleration (for reviews, see Blandford 
& Eichler 1987; Malkov & Drury 2001) provides a ba- 
sic framework for our understanding of cosmic-ray ac- 
celeration in coUisionlcss SNR shocks, though some key 
ingredients in the theory, especially the generation of 
magnetohydrodynamic (MHD) waves and their interac- 
tions with cosmic-ray particles (Bell & Lucek 2001; Bell 
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2004; Diamond & Malkov 2007), are yet to be understood. 

X-ray observations of shell-type SNRs can in principle 
probe the "microphysics" of the shock acceleration pro- 
cess; the development of MHD waves can be explored by 
measuring the energy of the spectral cutoff that is reg- 
ulated by the acceleration rate, which, in turn, is deter- 
mined by the shock velocity and the energy density of 
stochastic magnetic fields relative to ordered fields. In 
previous work, it has been presumed that the nonthermal 
X-ray spectrum is formed by a gradually steepening syn- 
chrotron spectrum extending from the radio band, as first 
outlined by Reynolds (1996). Experience shows that it is 
difficult to establish the presence of a cutoff (or roll off) 
in X-ray spectra that cover only a decade or so in energy, 
like data typically obtained by ASCA or Chandra. Secure 
detection of the expected spectral steepening requires a 
wider bandpass X-ray spectrum. In contrast, radio obser- 
vations of SNRs sample the energy distribution of acceler- 
ated electrons at considerably lower values (typically GeV 
energies) where the spectral form is more nearly power- 
law like, and where the power-law index in the framework 
of shock acceleration theory is closely related to the shock 
compression ratio. 

The SNR RX J1713. 7-3946 is a promising candidate 
to search for a spectral cutoff or steepening using mea- 
surements over as wide a wavelength coverage as pos- 
sible in the X-ray domain. RX J1713.7— 3946 (also 
known as G347.3— 0.5) exhibits the brightest nonther- 
mal X-rays among shell-type SNRs; the total nonther- 
mal X-ray flux (2-10 keV) integrated over the entire SNR 
amounts to ^ 5 x 10^^° erg cm~^ s~^, which is not ex- 
ceeded by any other shell-type SNR. This was the first 
SNR to be confirmed as a TeV 7-ray source using ground- 
based Cherenkov telescopes by two groups, CANGAROO 
(Muraishi et al. 2000; Enomoto et al. 2002) and H.E.S.S. 
(Aharonian et al. 2004; Aharonian et al. 2006; Aharonian 
et al. 2007). Detections of TeV 7-rays distributed over 
the entire remnant made with H.E.S.S. with a 1-10 TeV 
flux of 3.5 X 10~^^ erg cm~^ s~^, make this object well 
suited for exploring cosmic-ray acceleration in shell-type 
SNRs. The spectacular H.E.S.S. 7-ray image with spatial 
resolution of a few arc-minutes appears generally quite 
similar to the X-ray image, suggesting an intimate phys- 
ical connection between the nonthermal X-rays and TeV 
7-rays. 

Due to its faintness in the radio band, the remnant 
eluded detection in early radio surveys and was only dis- 
covered thanks to X-ray observations during the ROSAT 
All-Sky Survey (Pfeffermann & Aschenbach 1996). ASCA 
revealed intense synchrotron X-ray emission in the north- 
western (NW) part of the remnant (Koyama et al. 1997). 
The spectrum here is best-fitted by a power law of 
photon index T = 2.4 ±0.1 with interstellar absorption 
A^H = (0.81 ±0.06) X 10^^ cm^^, and is widely presumed 
to be synchrotron in nature. Diffuse synchrotron X- 
ray emission appears to be distributed over the entire 
face of the SNR with similar index (F ~ 2.3 ± 0.2) and 
without measurable thermal X-ray emission (Slane et al. 
1999). Faint synchrotron radio emission was detected 



in RX J1713.7— 3946, especially toward the western shell 
(Slane et al. 1999; Lazendic et al. 2004). With Chandra 
observations, it was found that the brightest part of the 
northwest X-ray shell consists of a complex network of 
bright filaments and knots embedded in more diffuse emis- 
sion (Uchiyama et al. 2003). The X-ray spectra of vari- 
ous components in the NW rim are similar to each other 
and well fit by a power-law function with F ~ 2.1-2.5 
(Uchiyama et al. 2003; Lazendic et al. 2004). Similar spa- 
tial and spectral characteristics were found in the south- 
western (SW) rim as well based on XMM-Newton data 
(Cassam-Chenai' et al. 2004; Hiraga et al. 2005). 

A nearby molecular cloud at ~ 1 kpc is likely associated 
with RX J1713. 7-3946, based on high sensitivity CO ob- 
servations made with the NANTEN telescope (Fukui et al. 
2003; Moriguchi et al. 2005). The association between the 
molecular cloud at ~ 1 kpc and the X-ray remnant is sug- 
gested by their close morphological match as well as by the 
detection of broad- line CO emission (Fukui et al. 2003). 
This is further strengthened by a detailed comparison be- 
tween the X-ray and molecular emissions (Uchiyama et 
al. 2005). The atomic and molecular hydrogen column 
density as a function of distance along the line-of-sight 
to the SNR reaches a value matching the X-ray absorb- 
ing column density at a distance of ~ 1 kpc (Koo et al. 
2004), thus providing independent support for an associ- 
ation with the nearby (1 kpc) molecular cloud. With this 
distance, the guest star that appeared in A.D. 393 possibly 
corresponds to the supernova explosion from which SNR 
RX J1713. 7-3946 originated (Wang et al. 1997). This 
association would make the remnant 1600 yr old. With 
this assumption, an X-ray point source (without an op- 
tical/radio counterpart) near the center of the remnant, 
IWGA J1713.4— 3949, is thought to be an associated neu- 
tron star (Slane et al. 1999; Lazendic et al. 2003), indi- 
cating a core-collapse supernova. It has been suggested 
(Slane et al. 1999; Ellison et al. 2001) that the supernova 
explosion occurred within the wind-blown bubble created 
by the massive progenitor star. We take the distance to 
SNR RX J1713. 7-3946 to be 1 kpc in this paper. 

Here we report on a wide band X-ray observation from 
0.4 to 40 keV of the SW part of RX J1713.7-3946 made 
with the Suzaku satellite. The high sensitivity in the 
12-50 keV band offered by the Hard X-ray Detector 
(HXD) on board Suzaku together with the X-ray Imaging 
Spectrometer (XIS) in the 0.4-12 keV band, allow us to 
measure for the first time the detailed spectral shape of 
the nonthermal X-ray emission from a supernova rem- 
nant. Since this is the first clear example of detecting ex- 
tended hard X-ray emission with the HXD, we performed 
an in-depth analysis of the hard X-ray measurements. In 
Section 2, we describe the observations and methods of ob- 
taining the data sets. In Section 3, we explain the analysis 
procedures to extract spectra and present spectral fitting 
for each of the HXD and XIS in turn. We then proceed 
with the implementation of HXD+XIS simultaneous mod- 
eling. We discuss the broadband X-ray spectrum in terms 
of the theory of diffusive shock acceleration in Section 4. 
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2. Observation and Data Reduction 

Wc observed the southwestern (SW) part of RX 
J1713. 7-3946 with Suzaku on 2005 September 26 (in 
the Performance Verification phase) for a total on-source 
duration of about 60 ks. The briglitest portion of the 
remnant is its northwest shell, which lies nearly on the 
Galactic plane, where we were concerned that an unex- 
pected transient source could appear and contaminate 
the HXD field of view. We therefore chose the second 
brightest part, i.e., the SW, to attempt a secure mea- 
surement of hard X-ray emission with the HXD. The 
X-ray Observatory Suzaku (Mitsuda et al. 2007), devel- 
oped jointly by Japan and the US, has a scientific pay- 
load consisting of two kinds of co-aligned instruments, 
the XIS (Koyama et al. 2007) and HXD (Takahashi et al. 
2007; Kokubun et al. 2007). The non-imaging HXD covers 
the hard X-ray bandpass of 10-600 keV with PIN silicon 
diodes (10-60 keV) and GSO scintillator (40-600 keV), 
both located inside an active BGO shield. We here make 
use of spectral data from the PIN diodes. The low back- 
ground of the HXD PIN enables us to study high energy 
emission from RX J1713.7— 3946 up to several tens of keV. 
Additionally, a narrow FOV of - 25' x 25' (FWHM) , which 
is defined by passive shields inserted in the BGO well-type 
collimator above the PIN diodes, is better suited to the 
study of extended emission in comparison with previous 
missions having larger FOVs. The XIS consists of four X- 
ray imaging CCD cameras covering a FOV of 17.'8 x 17'.8; 
three are front-illuminated (FI: 0.4-12 keV), one is back 
side illuminated (BI: 0.2-12 keV), and all have good en- 
ergy resolution of ~ 130 eV at 6 keV. 

The observation log for our observations is shown in 
Table 1. In order to estimate possible contributions from 
diffuse emission along the Milky Way in the hard 10- 
60 keV band, we observed two nearby regions in which 
no apparent hard X-ray point sources exist (off-pointings, 
OFFl and 0FF2). The three fields of view are depicted in 
Figure 1, overlaid on the ASCA 1-5 keV image (Uchiyama 
et al. 2005). More recently, one year after the PV phase 
observations discussed here, we have covered almost the 
entire remnant using 10 separate pointings during the 
Suzaku AO-1 period. These results will be published else- 
where (Tanaka et al., in preparation). 

We used data products made with version 1.2 of the 
pipeline processing. For the XIS analysis, we retrieved 
"cleaned event files" which were screened with standard 
event selections. We further screened the events with the 
following criteria: (1) cut-off rigidity larger than 8 GV 
and (2) elevation angle from the Earth rim larger than 10°. 
For HXD data, "uncleaned event files" were screened with 
standard event screening criteria. We excluded events 
during SAA passages and Earth occultation, and also 
those with cut-off rigidity less than 8 GV. The effective ex- 
posure times after these screenings are shown in Table 1. 
In this paper, data analysis was performed with HEAD AS 
6.0.1. 
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Fig. 1. Suzaku field of view for on-source (green) and off— 
source pointings (red), overlaid on the ASCA GIS image (1-5 
keV) taken from Uchiyama et al. (2005). North is up and 
east is to the left. Galactic coordinates are also shown. The 
inner square corresponds to the FOV of the XIS observations 
and the outer diamond shows the rough shape of the FOV 
(50% effective area) of the PIN detectors (see Fig. 6). The 
dashed square for the on-source demonstrates the full FOV of 
the PIN. 

3. Analysis 

S.l. HXD Data Analysis 

The HXD PIN achieves the lowest background level 
among all previous hard X-ray instruments (Kokubun et 
al. 2007). Also a rather narrow FOV defined by a passive 
collimator suppresses contamination from nearby discrete 
sources and diffuse emission such as the Galactic ridge 
emission (Valinia & Marshall 1998) and the Cosmic X- 
ray background (CXB). However, even for a bright object 
like this target, the PIN spectrum is still exceeded by the 
background — time-variable instrumental background in- 
duced in one way or another by cosmic-rays and trapped 
charged particles in orbit. The HXD instrument team 
has developed an effective method (Watanabe et al. 2007) 
of modeling the time-dependent non- X-ray background 
(NXB) by making use of the PIN upper discrimina- 
tor (UD) signal that monitors passing charged particles 
through the silicon PIN diode. The background spectrum 
is generated based on a database of NXB observations 
accumulated to date during night- and day-earth observa- 
tions. Note that in our analysis the background spectral 
model was generated with 100 times the actual number of 
background counts in order to minimize the photon noise 
on the background. The current NXB model is shown to 
be accurate within ~ 5% (Mizuno et al. 2006). In Figure 2, 
we show the light curve of the raw count rate in the 12-40 
keV range as a function of time for the on-source obser- 
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Table 1. Summary of the Suzaku observations of RX J1713.7— 3946 



Pointing 


Obs. ID 


Coord. (J2000) 
RA, DEC 


Exposure 
XIS/HXD 
[ks] 


Date 


RX J1713. 7-3946 SW 


100026010 


17'>12™17f0, -39°56'11" 


55/48 


26/9/2005 


OFFl 
0FF2 


100026020 
100026030 


17'^09'"31?9, -38°49'24" 
17'^09'"05?1, -41°02'07" 


28/24 
28/28 


25/9/2005 
28/9/2005 
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Fig. 2. HXD PIN lightcurves in the 12-40 keV energy band. 
The black, red, and green points show the raw, background 
(NXB + CXB), and background-subtracted data, respec- 
tively. The NXB curve is constructed with the background 
model of Watanabe et al. (2007). The background-subtracted 
count rate is almost constant within statistical errors, as ex- 
pected for the extended object like this, except for high back- 
ground periods. 

vation, the NXB count rate constructed with the model, 
and the background-subtracted lightcurve. It should be 
noted that the background-subtracted lightcurve has to 
be carefully examined in order to assure the validity 
of the background estimation. Based on visual inspec- 
tion, we discarded time regions when the total count rate 
was larger than 0.6 counts s^^ since the background- 
subtracted lightcurve shows an enhancement in those time 
regions. (When fitted with a power law, the difference of 
photon index for fits with and without these time regions 
is AF ~ 0.2.) The source intensity after background sub- 
traction is about 50% of the NXB level, implying that the 
systematic error on the hard X-ray flux due to a miss- 
estimation of the background should be about 10%. 

The HXD PIN spectra with and without background 
subtraction are shown in Figure 3 together with the back- 
ground model. Since the NXB model does not include 
contributions from the CXB, a simulated spectrum of the 
CXB was added to the NXB model. Following Gruber 
et al. (1999), who reanalyzed HEAO-1 observations in 
the 1970's, our spectral model for the extragalactic back- 
ground radiation was chosen as 

dN _ „ _i 9Q / CkcV 
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Fig. 3. Hard X-ray spectrum of the southwest part of 
RX J1713.7-3946 recorded with the HXD PIN. The data 
points in black show the raw spectrum, red points represent 
the background model (NXB+CXB), and green points are the 
background-subtracted spectrum. 

where ekcV = e/1 keV and Cp = 41.13. We estimated the 
expected CXB signal in the HXD-PIN spectrum using 
the latest response matrix for spatially uniform emission, 
ae_hxd_pinf lat_20060809 . rsp. The contribution from 
the CXB flux is estimated to be ~ 5% of the NXB, com- 
parable to the current systematic uncertainty of the NXB 
model itself. We confidently detect hard X-rays from 
RX J1713. 7-3946 up to 40 keV; at higher energies back- 
ground uncertainty dominates the source flux. 

The PIN data from the two OFF-pointing observations 
were analyzed in the same way as those from the on- 
source pointing. After subtraction of NXB and CXB, 
we do not find any significant hard X-ray emission that 
exceeds the background uncertainty. If we fit an NXB- 
subtracted spectrum of the OFF pointings (OFFl and 2 
are added) with the CXB spectral form defined by Eq. (1) 
but varying its normalization, we find a normalization fac- 
tor of 1.7 ±0.2 (where the uncertainty includes only sta- 
tistical errors). This suggests the possible presence of dif- 
fuse emission from the Galactic plane (Valinia & Marshall 
1998) at a level of about 70% of the CXB. This is below the 
uncertainty in the NXB model and will not significantly 
affect the RX J1713.7— 3946 spectrum. It should be noted 
that the RXTE PGA spectrum of RX J1713. 7-3946 in 
^)the 2-30 keV range (Pannuti et al. 2003) is likely contam- 
inated by the Galactic ridge emission because of its much 
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Fig. 4. Background-subtracted PIN spectra, using the back- 
ground model (black) and off-source spectrum (red). 

larger FOV 1° FWHM). Instead of the background 
model, we subtracted the OFF-pomting spectrum from 
the on-source spectrum as shown in Figure 4. We found 
that the two methods of background subtraction gave al- 
most identical results within the statistical error, though 
the OFF-subtracted spectrum seems to be systematically 
lower than the model-subtracted spectrum. This could in- 
dicate the level of the systematic error in the background 
model or some contribution from the diffuse ridge compo- 
nent. 

Before proceeding with the spectral fitting, we checked 
the detector response against an extended source like 
RX J1713.7— 3946, since this is the first clear case that 
the Suzaku HXD gives a spectral measurement for any 
source larger than the PIN collimator. We performed 
Monte Carlo simulations with the code simHXD (Terada 
et al. 2005) to study the PIN response for an extended 
source. Hard X-rays from a 2° x 2° sky area centered on 
the optical axis were uniformly distributed over the en- 
trance aperture of the HXD, and their interactions with 
the detector system were tracked. Figure 5 shows the 
simulated spectra of a diffuse source (red) and a point 
source at the XIS- nominal position (black). The com- 
parison demonstrates that spatial extent does not affect 
the observed spectral shape, consistent with expectations 
that the built-in passive collimator (50 /Ltm thick phosphor 
bronze sheets) should almost completely absorb hard X- 
rays below 80 keV and, consequently, that the angular re- 
sponse is not energy-dependent. In Figure 6 we show the 
angular response of the PIN detector determined by the 
passive collimator. The difference of the spectral shape is 
so small that it corresponds to the difference of photon in- 
dex of 0.01 when fitted with a power-law model using the 
same response matrix. We conclude that the extended na- 
ture of the source does not cause any spectral steepening 
or flattening, at least for energies below 40 keV. In the 
following, we therefore fit the PIN spectrum simply us- 
ing the point-source response matrix at the XlS-nominal 
position, aejixd_pinxinom_20060814 . rsp. 
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Fig. 5. Simulated PIN spectra for a point (red) and extended 
source (black). The bottom panel shows the ratio between 
them. This demonstrates that the energy-dependence of the 
PIN response for an extended source is same as that for a 
point source. 
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Fig. 6. Angular response of the PIN detector as defined by 
the built-in passive collimator. Contours indicate the effective 
area as a function of incoming photon direction. The angular 
response is independent of energy below 80 keV. 
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We fitted the background-subtracted PIN spectrum 
with a simple power law model: dN/de oc e^'". The best- 
fit parameters are summarized in Table 2. We present 
two cases, one using the NXB-f-CXB model as a back- 
ground spectrum, and the other using the OFF-pointings 
as a background template. Both methods agree with 
each other. Hereafter we only use spectral results with 
the NXB+CXB model. The photon index obtained, F = 
3.2 ±0.2, is much larger than that for previous soft X-ray 
(below 10 keV) results, F~ 2.4 (Slane et al. 1999). Given 
the current reproducibility of the NXB model of ~ 5%, the 
systematic error due to uncertainties in the NXB modeling 
is found to be AF ~ 0.2. A spectral cutoff around 10 keV 
seems necessary to account for the spectral transition from 
the flat X-ray spectrum with F = 2.4 ±0.05 (Slane et al. 
1999) into the steep spectrum with F = 3.2 ± 0.2 in the 
PIN domain. Below we further confirm this important 
finding based on spectroscopy with the XIS. 

3.2. XIS data Analysis 

Figure 7 shows the XIS images from our observation in 
the soft (1-5 keV) and hard (5-10 keV) bands. Both im- 
ages were smoothed with a Gaussian kernel with cr = 0'.3 
and background-subtraction was carried out for each of 
them. We utilized the Night Earth Background Database 
consisting of event data accumulated while the satellite 
was viewing the night earth where the non X-ray back- 
ground becomes dominant. After background subtraction, 
the vignetting effects of the X-ray mirrors were corrected 
by means of the simulation program xissim (Ishisaki et 
al. 2007). It should be noted that the two XIS images in 
the soft and hard bands are similar to each other. This in- 
dicates that spectral changes reported by Cassam-Chenai' 
et al. (2004) across the SW part are not dramatic in our 
data set. 

The XIS provides higher quality spectra of the SNR 
RX J1713.7— 3946 compared with previous missions 
thanks to its low background level, good energy resolu- 
tion, excellent energy response below 1 keV, and large 
effective area (Koyama et al. 2007). We accumulated 
X-ray photons from the region shown in Fig. 7. Since 
RX J1713.7— 3946 is situated near the Galactic plane, the 
background spectrum needs to be accumulated from a 
nearby, source-free region. We extracted the background 
spectra from the OFF observations. In the fitting pro- 
cedure, the standard RMF files (response matrices of the 
XIS) version 2006-02-13 were used, whereas ARF files (re- 
sponse of the XRT) were produced with the xissimarf gen 
software version 2006-10-26 (Ishisaki et al. 2007). The en- 
ergy range of 0.4-12 keV, to which the XIS FI chips are 
sensitive, is used in fitting, but a range of 1.7-1.9 keV is ex- 
cluded from the spectra because of large systematic uncer- 
tainties in the response matrices in this band. Background 
uncertainty is most significant for the low energy (<0.6 
keV) part of the spectrum where the signal to background 
ratio drops. We do not take account of this issue in detail, 
except to note that the spectral fits reported here are not 
affected by this. Finally, we co-added the three spectra, 
RMF files, and ARF files from the FI chips to produce a 



single data set for the XIS. 

Following previous studies, we first fitted the XIS spec- 
trum in the 0.4-12 kcV band with a simple power law 
attenuated by interstellar absorption. Wc found a power 
law model fails to fit the data, giving a large reduced 
chi-squared of xl = 1-81 for 302 degrees of freedom. 
Though statistically unacceptable, the derived photon in- 
dex of F = 2. 39 ±0.01 and absorbing column density of 
~ (0.87 ±0.01) X 10^^ cm~^ are in good agreement 
with the ASCA result in this region: F 2.40 ± 0.05 and 
Nu = (0.79 ±0.50) X 10^2 cm'^ (Slane et al. 1999). In 
Figure 8 (left panels) we show the XIS spectrum with 
the best-fit power-law function and residuals (data mi- 
nus model). There remains a correlated pattern to the 
residuals in the power-law fit. This failure of the power- 
law fit argues that a more complicated spectral model is 
required. Furthermore, the difference in fitted photon in- 
dices between the XIS and HXD PIN data says that the 
model needs to roll-off or steepen with energy across the 
broad bandpass. 

We therefore fit the XIS spectrum with three other 
models: a power-law with two types of an exponential 
cutoff (cutoff power law), (1) exp [— (e/ec)] and (2) 

e~'"exp —\ftJTc , and (3) a broken power law. Unlike 

the pure power-law case, these mathematical functions 
give acceptable fits (see Table 3). Figure 8 (right panels) 
shows the XIS spectrum with the best-fit e~'"exp [— (e/ec)] 
function together with its residuals. The wavy residuals 
seen in the power-law fit clearly disappear by introduc- 
ing the cutoff. The cutoff power-law fit gives a power law 
F = 1.96 ± 0.05 with a cutoff at = 8.9tJ:g keV, where 
the fit function has a form e^'"exp[— (e/ec)]. Also, we find 
a hard power law F = 1.50 ± 0.09 with a lower cutoff at 



. = 1.2+n-2 keV with the fit function of e 



The broken power-law model also reduces the correlated 
pattern in the residuals. The best-fit function is com- 
posed of two distinct power laws of Fi = 2.18 ± 0.03 and 
F2 = 2.53 ± 0.02 with a break energy of = 3.13 ± 0.15 
keV. 

We also performed spectral fits with the so-called 
SRCUT model (Reynolds 1998) in the XSPEC package, 
which describes the synchrotron spectrum from electrons 
following an exponentially-cutoff power law distribution 
in energy. In XSPEC, its spectral shapes for different pa- 
rameters are tabulated based on numerical calculations 
convoluting a well-known synchrotron formula for a single 
electron with the energy distribution of electrons of the 
form E^''^ cxp{~E/ Eq). The synchrotron spectrum has 
a power-law form at radio frequencies (with photon index 
F = a + 1 , where a is nominally the spectral index in the ra- 
dio band), with a rolloff energy eioii (nominally the energy 
where the exponential cut-off reduces the fiux by a factor 
of 4 below the extrapolation of the radio power-law spec- 
trum). Wc fixed the index to a typical value for SNRs of 
a = 0.5. A good fit was obtained with a best-fit rolloff en- 
ergy of eroii = 0.95 ±0.04 keV. The roll-off energy obtained 
in this way is similar to the cutoff energy Ec = 1-2 keV ob- 
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Broad-band X-ray spectrum of SNR RX J1713. 7-3946 

Table 2. Power-law fitting to the HXD PIN spectrum* 



Pointing 


Background 


r 


Flux T 










(10-40 keV) 




RX J1713. 7-3946 SW 


model (NXB-hCXB) 


3.2±0.2 


2.5±0.1 


1.15 (36) 


RX J1713. 7-3946 SW 


0FFl-h2 


3.2±0.3 


2.3±0.2 


0.92 (27) 



* Errors represent 90% confidence, 
t In units of mCrab. 





Fig. 7. The XIS images of RX J1713. 7-3946 in (a) the soft 1-5 keV and (b) the hard 5-10 keV bands. Three FI chips are combined. 
North is up and east is to the left. The images have dimension 17' X 17'. The color scale indicates the surface brightness in a linear 
scale. The green box shown in (a) is the region from which we extracted the XIS spectrum. 



tained for a spectral form e exp 



This is rea- 



sonable because the SRCUT model has an approximate 
form e~^exp —-^J e/croii in a narrow bandpass around 
the cutoff region. 

3.3. Suzaku Wide Band Spectrum 

We now combine the HXD PIN spectrum with the 
XIS spectrum to fully exploit the broadband capability of 
Suzaku. Since the source is extended and the HXD PIN 
observes a much wider field than the XIS, the relative nor- 
malization factor derived for a point source, k= 1.13±0.01 
(Ishida et al. 2006), is not applicable in our case. We note 
that the relative normalization factor k is defined here as 
the flux density at 12 keV determined by the PIN relative 
to that determined by the XIS. Therefore, it is necessary 
to calculate the proper scaling factor between the XIS and 
the HXD for a diffuse emission source. We determined the 
normalization factor k, which takes into account the larger 
FOV of the PIN, by convolving the emission distribution 
of the ASCA image with the angular response of the HXD- 
PIN defined by the passive collimator. It is assumed that 
the spatial distribution of X-ray surface brightness is in- 
dependent of energy in the ASCA-Suzaku band over the 
entire remnant. We determined the normalization factor 



to be K = 6.6 using the ASCA 5-10 keV image, with a 
possible systematic error of ~ 10-20%. (The value of k 
does not depend sensitively on the choice of ASCA band- 
pass. If instead we use the ASCA 1-5 kcV image a slightly 
different value, k = 6.5, is obtained.) 

In Figure 9, we show the XIS-I-PIN spectra together 
with the power law, cutoff power laws, broken power law, 
and SRCUT models. The model parameters are the best- 
fit values determined by fitting the XIS spectrum alone 
(see Table 3), and the best- fit models are extrapolated 
to the PIN bandpass after multiplication by the factor 
K ~ 6.6. The power-law model clearly conflicts with the 
PIN data as does the broken power law model, which 
was a good description of the XIS data alone. Although 
the shape of the SRCUT model in the PIN band (which 
can be characterized as a power law with a photon in- 
dex of r ~ 3.3) is a good match to the data, the level 
of hard X-rays this model produces is about 40% too 
high, somewhat more than the estimated systematic er- 
ror in the normalization factor k of ~ 20%. If we change 
the exponent, a, of the SRCUT model by ±0.1, this re- 
sult still holds. Although the cutoff power-law function 
of e^'" exp [— (e/ec)] satisfactorily reproduces the level of 
the observed PIN flux, the model flux in this case (ap- 
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Fig. 8. (Left) The XIS spectrum in the 0.4-12 keV band together with the best-fit power-law model. The lower panel shows the 
residuals between the data and model. (Right) Same as the left panel but with the best-fit cutoff power-law model. 

Table 3. Model fitting to the XIS spectrum (0.4-12 keV)* 



Function 



power law (e ^ ) 



(10^ 



(keV) 



0.87±0.01 2.39±0.01 

0.77 ±0.01 1.96 ±0.05 

0.74 ±0.02 1.50 ±0.09 

0.79±0.01 2.18 ±0.03/2.53 ±0.02 3.13±0.15 



(10- 



Fl-io keV 



e ^ exp 
exp 



1 r, + 0.3 



broken power law 
srcut 



0.77 ±0.01 



1.5 fixed 



Notes. 

* Errors quoted at 90% confidence, 
t Corrected for absorption. 



0.95 ±0.04 



6.07 ± 0.09 

'-'•"'-0.05 

5.61 ±0.40 

5.75±0.17 
5.71 ±0.02 



1.81 (302) 
1.03 (301) 

1.00 (301) 

1.00 (300) 
1.05 (302) 



proximately a power law with T i 
falls more steeply than the data, 
form of e^'^exp 



3.7 in the PIN band) 
A slower exponential 
, which gives F ~ 3.5 in the PIN 

band for its best-fit model, seems to yield the best match 
to the PIN data both in the absolute flux and spectral 
shape. We conclude that the evidence for a curvature in 
the broad band Suzaku spectrum of RX J1713.7— 3946 is 
secure. However, it is premature at this time to argue 
that, for example, the cutoff power-law models are a bet- 
ter description of the emission conditions in the SNR than 
the SRCUT model based on the XIS+PIN joint spectral 
analysis, simply because the two data sets do not come 
from the same spatial regions. Modest spatial variations 
in the emitting conditions could produce flux or spec- 
tral differences similar to those seen in the PIN band for 
the cutoff power-law and SRCUT models. Furthermore, 
even spatial/spectral variations in the synchrotron emis- 
sion within the joint XIS±PIN field of view (as seen by 
Cassam-Chenai et al. (2004)) could account for the mod- 
est differences in the cutoff power-law and SRCUT models 
in the PIN band. 

3.4- Upper Limits on Thermal Emission 

One of the most remarkable characteristics of 
RX J1713.7— 3946 is the complete lack of a thermal X-ray 
emission component as expected from the hot shocked gas 
in the SNR. Our XIS spectrum, similar to previous stud- 



ies (Slanc et al. 1999; Cassam-Chenai et al. 2004), shows 
no signs of thermal emission. In the following we set an 
upper limit on the density of shocked hot gas from the 
upper limit on the thermal emission measure. 

To derive the allowable level of thermal emission un- 
derneath the synchrotron spectrum, we added a ther- 
mal component to the cutoff power law model obtained 
by fitting the XIS spectrum (see Table 3 for the pa- 
rameters). We used the APEC model for the thermal 
emission model and assumed solar composition. Then, 
we determined 3(7 upper limits on the emission measure, 
EM = / rientidV / AttD"^ , as a function of electron temper- 
ature kTc over the range 0.05 keV to 1 keV. Above 1 keV, 
the upper limit has only a weak dependence of tempera- 
ture (Slanc et al. 1999). Here, ric, tt-h, V , and D are the 
number density of electrons and hydrogen, the volume of 
the hot gas, and the distance to the remnant, respectively. 

Figure 10 shows the Za upper limits as a function 
of temperature. The upper limits grow quite large for 
fcTc ^ 0.2 keV. The emission measure can be as large 
as 10"'^^ cm~^ if we assume a temperature as low as 
~ 0.05 keV. From the emission measure, we can esti- 
mate the number density of shocked gas. By assuming 
Uc = 1.2riH, uniform distribution of gas inside the rem- 
nant, and that the spectral region covers 8% of the whole 
volume of the remnant (a sphere of 30' radius), we relate 
the gas density (71 = tih) to the emission measure as 



No. ] Broad-band X-ray spectrum of SNR RX J1713. 7-3946 




S d 
3 to 

o a 



1 ' ' 1 • . . 1 1 1 

: (d) -v.. Broken P 

^ / ^\ 


L 1 

^ 


|i ji |i 1 1 1 1 1 — 1 — 1 — 1 ' ' ' 1 1 

1 . 1 





1 10 
chaimel energy (keV) 



channel energy (keV) 




1 10 
channel energy (keV) 




chaimel energy (keV) 




1 10 
channel energy (keV) 



Fig. 9. XIS-I-PIN spectra in the energy range 0.4—40 keV are shown with the best-fit models obtained by fitting the XIS spectrum 



alone: (a) power law, (b) e exp [— (e/tc)] (c) e '"exp 
the ratio between the date points and the model values. 



(d) broken power-law, and (e) SRCUT. The lower panels show 
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Fig. 10. The Scr upper limit on thermal emission 
in RX J1713.7— 3946 expressed as emission measure 
(EM = J nenndV/i-n-D^) vs. electron temperature fcTg. 



71=1 



EM 



101 



1/2 



D 



1 kpc 



-1/2 



-3 



(2) 



For a distance of 13 = 1 kpc. the upper limit on the gas den- 



sity ranges from 9 x 10 



to 2 cm when we vary 



the temperature from 1 keV to 0.05 keV. This is broadly 
consistent with previous estimates based on ASCA (Slane 
et al. 1999) and XMM-Newton (Cassam-Chenai et al. 
2004). 

4. Discussion and Conclusions 

As established in previous studies (Koyama et al. 
1997; Slane et al. 1999; Uchiyama et al. 2003; Lazendic et 
al. 2004; Cassam-Chenai et al. 2004; Hfraga et al. 2005), 
the X-ray spectrum extracted fr'om any portion of RX 
J1713.7— 3946 shows a featureless continuum, for which 
an absorbed power law offers a reasonable fit. The ASCA 
spectrum of the SW part was characterized by a photon 
index T = 2.40 ± 0.05 (Slane et al. 1999). 

The Suzaku broadband spectrum for the SW part of 
RX J1713.7— 3946 reported here clearly indicates the pres- 
ence of a spectral cutoff as it shows a spectral transition 
from a relatively flat power law (F ~ 2.3) in the soft X- 
rays to a steep power law (F ^ 3.3) at harder photon ener- 
gies. In fact, the nonthermal featureless X-ray spectrum of 
RX J1713.7— 3946, as well as other nonthermal remnants 
such as SN 1006 (Koyama et al. 1995) and Vela Jr. (Slane 
et al. 2001), are widely thought to arise from the gradually 
steepening high-energy part of the synchrotron radiation 
spectrum produced by TeV-scale electrons, which are ac- 
celerated through the mechanism of diffusive shock accel- 
eration (Reynolds 1998). The cutoff in the Suzaku spec- 
trum provides firm observational support for this widely- 
accepted presumption. We found that the cutoff energy, 
presumably set by synchrotron cooling, appears differ- 
ent from one spectral model to another, in a range of 
Ecutofi ^ 1^9 keV. The higher cutoff, ec ~ 9 keV, is based on 



the power law with exp(— e)-type cutoff model, while the 
lower bound is derived with exp(— ^/e)-type cutoff. Given 
this fact, wc need to be careful in discussing the physical 
meanings of the position of the spectral cutoff. Wc briefly 
explore the implications of the Suzaku broadband spec- 
trum in a separate paper (Uchiyama et al. 2007). More 
detailed discussions on this issue will be given elsewhere 
(T. Tanaka et al. in preparation). 

If the magnetic field exceeds 10 /iG in the shell (e.g., see 
Berezhko, & Volk (2006)), the TeV 7-rays detected from 
the remnant should be explained by protons producing 
7r°-decay 7-rays, rather than by electrons emitting 7-rays 
via inverse Compton scattering (Aharonian et al. 2006). 
The hadronic model requires the total energy content of 
protons to amount to 



Wp«10 



50 



D 



1 kpc 



n 



1 cm 



erg. 



(3) 



in order to explain the observed TeV 7-ray flux 
(Aharonian et al. 2006). Therefore, a matter density of 
~ 0.2 cm~'^ is needed, assuming the typical kinetic en- 
ergy released by a supernova of 10^^ erg and a conversion 
efficiency to high energy protons of 50%. To reconcile 
this density value with the upper limit on thermal X-ray 
emission (§3.4), the electron temperature must be of or- 
der 0.1 keV (or lower), otherwise thermal emission (oc n^) 
would have been detected with the Suzaku XIS. The low- 
temperature of thermal electrons may indicate that effi- 
cient particle acceleration drains the energy of the shock, 
thus suppressing gas heating at the shock (Ellison et al. 
2001). 

Strong (high Mach number) shocks in young shell-type 
SNRs do accelerate high-energy particles with remarkably 
high efficiency, thus providing an ideal laboratory to study 
high-energy particle acceleration in the Universe. Suzaku 
has excellent capability for determining the spectral form 
of X-ray synchrotron emission thanks to its broadband 
spectral coverage from 0.4 to 40 keV. The observed cutoff 
energy in the synchrotron X-ray spectrum suggests that 
shock acceleration proceeds at nearly the maximum pos- 
sible rate. Combined with the broad TeV 7-ray spectrum 
(in the 0.3-100 TeV band) measured recently with HESS, 
our Suzaku observations of RX J1713.7— 3946 provide 
stringent constraints on shock-acceleration of the highest 
energy particles in our Galaxy. 
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